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Summary. — We survey theoretical and experimental/observational results on 
general-relativistic spin (rotation) effects in binary systems. A detailed discussion is 
given of the two-body Kepler problem and its first post-Newtonian generalization, 
including spin effects. Spin effects result from gravitational spin-orbit and spin-spin 
interactions (analogous to the corresponding case in quantum electrodynamics) and 
these effects are shown to manifest themselves in two ways: (a) precession of the spin- 
ning bodies per se and (b) precession of the orbit (which is further broke down into 
precessions of the argument of the periastron, the longitude of the ascending node 
and the inclination of the orbit). We also note the ambiguity that arises from use of 
the terminology frame-dragging, de Sitter precession and Lense-Thirring precession, 
in contrast to the unambiguous reference to spin-orbit and spin-spin precessions. 
Turning to one-body experiments, we discuss the recent results of the GP-B experi- 
ment, the Ciufolini-Pavlis Lageos experiment and lunar-laser ranging measurements 
(which actually involve three bodies). Two-body systems inevitably involve astro- 
nomical observations and we survey results obtained from the first binary pulsar 
system, a more recently discovered binary system and, finally, the highly significant 
discovery of a double-pulsar binary system. 
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2 R. F. O'CONNELL 

1. Introduction 

This manuscript is, in essence, a progress report on general-relativistic spin effects 
in binary systems, which I lectured on in Varenna summer schools in 1974 [T] and 1975 
[2]. Here, we use "spin" in the generic sense of meaning "internal spin" in the case 
of an elementary particle and "rotation" in the case of a macroscopic body. My 1974 
lecture concentrated on the one-body problem [3] and analysis of the Schiff gyroscope 
experiment [4J, now referred to as the Gravity Probe B (GP-B) experiment [5J, which 
was launched on April 20, 2004 [6]. In addition, I discussed the possibility of C, P, 
and T effects occuring in the gravitational interaction. Not long after the 1974 summer 
school, a momentous event occured: the discovery of a pulsar in a binary system by 
Hulse and Taylor [7] which bought into play a new array of possibilities for observing 
general relativistic effects. In particular, this motivated Barker and me to calculate 
two-body spin precession effects [HI H] , more details which were reported on later at the 
1975 Varenna summer school. Thus, other than summarizing the key results from these 
lectures, I will concentrate here on discussing progress since then, in theory, experiment 
and observation. Section 2 is devoted to a general discussion of the Kepler problem and 
its first post-Newtonian generalization (order c~^) including spin effects. In section 3, 
we discuss a one-body system and related experiments and observations, not only the 
GP-B experiment but also the more recent results of Ciufolini and Pavlis [10] and the 
very early theoretical work of de Sitter [TTJ [T^] and Lense-Thirring TJl [TJ as well as the 
related work of Williams et al. [15] and Murphy et al. [16: . Section 4 concentrates on 
2-body systems and we discuss, in particular, the observations of Stairs et al. [T7] and 
the momentous discovery of a double binary system [18j . 

2. The Kepler Problem and its Post-Newtonian Generalization 

For a binary system, the Newtonian gravitational potential gives rise, for energy 
E < 0, to elliptic motion. We are interested in post-Newtonian generalizations (or- 
der c~^) arising from both spin and general relativistic effects. First, we present our 
notation [5]. Let ttt-i, ri, vi. Pi, S'-^-', n^^^ J*^^) and a;^^) denote the mass, position, veloc- 
ity, momentum, spin, unit vector in the spin direction, moment of inertia and angular 
velocity, respectively, of body 1. The same symbols, with 1^2, denote the corre- 
sponding quantities for body 2. In the center-of-mass system we have r = ri — r2 and 
P = Pi = — P2 = /iv where v = Vi — V2. The reduced mass and the total mass are 
given by 

TTllTTlQ 

(2.1) fj.= , M = mi -1-7712. 

7711 + "T-2 

Also, the orbital angular momentum is 

(2.2) L = r X P = in. 
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where n is a unit vector in the L-dircction. In additon, 



(2.3) 



a2 



Lin _ /GMy _2tt _ 

(i-e2)3 " y^J ^Y^' 



where T is the orbital period, e is the eccentricity, a is the semi-major axis, and ui is the 
average orbital angular velocity. 

For a Newtonian elliptic orbit for two spherically symmetric bodies, the energy E, the 
orbital angular momentum L, and the Runge-Lcnz vector A are constants of the motion, 
which can be written as 



(2.4) E/^l - i«2 - GM/r, 



(2.5) L/fi = r X V, 

and 



(2.6) A/m = V X (r X v) - GMr/r = -(v x L) - GMr/r, 

M 

where the direction of A is along the major axis from the focus to the perihelion and is 
thus perpindicular to L. However, when spin and post-Newtonian effects are taken into 
account (order c"^ beyond Newtonian theory) we found that the secular results for the 
precession of the orbit are 



(2.7) Eav = 0, 



(2.8) 



(2.9) 
where 

(2.10) n* = n*(^) + n*'-^^ + fi*(2) + fi*(i'2)^ 
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consists of contributions from a spin-independent term ft*^^' (Einstein's precession of 
the perihelion), in addition to S^^^ S*^^^ and S^^^S*^^) terms, respectively, which are all 
given explicitly in [8]. We note that, on the average, L and A precess at the same rates, 
so that the Kepler ellipse as a whole precesses with the angular velocity il*. It was also 
found convenient (in order to adhere to astronomical and space physics practise) to write 
rj* in the form 



„^ dfl du) di nn X n 

(2.11) n=— no + — n' 



dt dt dt |no X n| ' 

where fl,ui and i denote the longitude of the ascending node, the argument of the peri- 
astron and the inclination of the orbit, respectively, in the reference system of the plane 
of the sky (the tangent plane to the celestial sphere at the center of mass of the binary 
system) [21 13]. In addition, no is a unit vector normal to the plane of the sky directed 
from the center of mass of the binary system towards the Earth. The angle between 
no and n is the inclination i. In the absence of spin, only the periastion precession is 
present. 

In addition, each spin precesses at rates 17^^-' and ^^^\ which when averaged over an 
orbital period we write as 



(2.12) SW^^WxSW, 

and similarly for 1-^2. However, as we recently stressed "gravitational effects due to 
rotation (spin) are best described, using the language of QED, as spin-orbit and spin- 
spin effects since they also denote the interactions by which such effects are measured; 
in fact these are the only such spin contributions to the basic Hamiltonian describing 
the gravitational two-body system with arbitrary masses, spins and quadruple moments. 
They manifest themselves in just two ways, spin and orbital precessions, and whereas 
these can be measured in a variety of ways (for example, as discussed above, orbital 
precession can be subdivided into periastron, nodal and inclination precessions) such 
different measurements are simply 'variations on the theme'" [19]. If fact, in the basic 
Hamiltonian for the problem [3, the spin-orbit terms are of the form S'-^-' • L and 8*^^^ • L 
whereas the spin-spin terms are of the form {[3(S*-^-' ■ r)(S^^^ • r)/r^ — S^^^ ■ S^^']}, which 
is completely analogous to what occurs in QED, as wc pointed out in [IJ. Thus, we were 
able to conclude [1 that all such effects in QED have their analogy in general relativity 
and, except for (important) numerical factors, the latter results may be obtained from 
the former by simply letting e^ -^ Gmim2. Hence it is convenient to write (with a similar 
equation for 1-^2) 



f2« = f^(i) -I- f2« 

av so ' ss 

(2.13) ^ n« + n^ll, 
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where the second equahty is often written to make contact with common usage, ds 
denoting de Sitter [TT1[12] and LT denoting Lense-Thirring [131 [14]. The former is also 
referred to as the geodetic effect and the latter as a "frame-dragging" effect. However, as 
we pointed out in [TO] , the de Sitter precession can also be regarded as a frame-dragging 
effect, which is even more obvious in the case of two-body systems (but which is also 
readily seen in the case of the GP-B experiment if one considers an observer in the frame 
of the gyroscope). Further confusion associated with the terminology "Lense-Thirring 
effect" arises from the fact that it is also used to discuss earth spin effects on orbital 
motion, as we discuss in more length below (subsect. 3.3). Thus, to avoid confusion, we 
feel that it is best to just use the terminolgy given in the first equality of (2.13). We also 
note that flyj and fl^^J (explicit results for which are given in [S]) do not depend on the 
spins and always point along L whereas n].J depends on the spin of body 2 and points 
along{[n(2)-3(n-n(2))n]}. 

To conclude this section, we make some general remarks. First, we note that since 
gravitational radiation only arises at order c~^, and defining the total angular momentum 



(2.14) J = L + S(i)+S(2\ 

we obtained [5] 



dJ 
(2.15) -^0. 

again underlining the inter-relationship between the spin and orbital precessions. Sec- 
ondly, we note that the largest precession is always the periastion precession and it is the 
only one that remains in the absence of spin. For a particular system, the magnitude of 
the periastion precession and the spin-orbit precessions depends significantly on the value 
of the product of the orbital angular velocity Co and the gravitational coupling constant 



(2.16) ag = (GM/c^a). 

where a is the semi-major axis of the Keplerian ellipse. Thus, for the GP-B satellite 
orbiting the earth, Q = 1.074 x 10"^s~i and ag = 6.958 x 10"^°, so that agCJ = 7.473 x 
lO^^"^.?"^ = (l.35 X 10^'*) /yr. On the other hand, for the double binary system, ui ~ 
7.112 X IQ-^s-^ and Ug = 4.348 x lO'^ so that agCu = 3.092 x lO^^s^^ = 5°.591/j/r. As 
we will see below, this difference, amounting to 4 orders of magnitude, is reflected in the 
significantly larger precessions (both orbital and spin) obtained for binary star systems 
compared to earth related systems. For example, this order-of-magnitude calculation is 
borne out by an exact calculation of the spin-orbit precession, which gives (in degrees 
per year) 1.84 x 10"'^ for the GP-B gyroscope and 4.8 and 5.1 for pulsars A and B, 
respectively, which constitute the double binary system, to be discussed below. Thirdly, 
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the inter-relationship between the spin and orbital precessions, as expressed in (j2.15p . 
may be seen explicitly by writing our exact basic results [8J in the following succient form 
(recalling again that the precession rates are averaged over an orbital period) 



dt 



(2.17) 



3 _ TO2+^/3 / (1) 



2(1 -e2) 






(n(^ 



(2)xS(i) -3(n.n(2)UnxS(i) 



dS(2) 
dt 



= te)+o(2)Us(2 



(2.18) 
and 



so ' ss 



TOi + fi/3 



2(1 - e2) "^^ M 

1 ^(1) 



2(1 - e2) ' Ma 



nxS(2)) 



nWxS(2))-3(n.nW)(nxS(2) 



'dt 



(2.19) 



dt 



dt 



-{(0«xSW 
'dS(i) dS(2 



J6„„ X 0„„ 



0« X S(i; 



J7 



(2) 



S(2)^ 



)]} 



which is consistent with (|2.14p and (|2.15p . This is another manifestation of the fact that 
there are simply two basic gravito-magnetic interactions in general relativity, spin-orbit 
and spin-spin. We note that fl*^ ' does not contribute to ^ whereas for ^ it is the 
dominant Einstein-Robertson periastion precession term. It is also immediately clear 
from these results that the various precession rates are proportional to ag times the 
average orbital rate, consistent with the fact that we are working to order c~^ (first post- 
Newtonian approximation). After this general discourse, we now turn to a discussion of 
various experiments to detect spin-orbit and spin-spin effects in gravitational theory. 



3. One-Body System 

3'1. The GP-B experiment. - As already mentioned, it was launched on April 20, 
2004. Data collection covered the period Aug. 2004-Sept. 2005, when the He in the 
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dewar had completely leaked out but, due to unexpected "-measurement uncertainty-," 
the announcement of final results has been delayed until December 2007. 

Here m2 (earth mass) >> nii (gyroscope mass) and S*-^-* is along the earth's rotation 
axis. Since fl^^^ is always along L, it turns out that a polar orbit results in flyj and 
illJ being at right angles resulting, in principle, in a clear separation of the two effects. 
However, due to a variety of other effects pnil21j . the analysis turns out to be much more 
complicated. The calculated precessional rates for the GP-B experiment, in units of arc- 
seconds per year (2.78 x 10"''' degrees per year) are 6.61 for fiio and 4.1 x 10~^ for Qls 
and the hope is that the accuracy achieved is better than 5 x 10""* (1.4 x 10~^ degrees) 
per year [22] . However, whether this is achievable in view of the problems encountered 
[HI nulls] remains to be seen. Among the concerns is the polhode motion arising from the 
fact that the gyroscopes are not perfect sphere. Thus, not only can a quadrupole moment 
arise from the manufacturing process but, in addition, even an ideal perfectly spherical 
gyro acquires a quadrupole moment due to its rotation [5D]. Using the quoted values of 
the gyro diameter and average spin rate [22] we obtain a value for (AI/I), arising from this 
effect, of 7.1 x 10~^, which is comparable to that arising from the manufacturing process. 
An additional problem as far as the experiment is concerned is that the two different 
quadrupole moments are located at an unknown angle from each other, aggravated by 
the fact that "-there was no way to know about which axis each gyro (of the four total) 
was initially spun up-" [23] . Moreover, the "-wide variation in spin-down rates-" will 
be reflected by a time variation in the value of (AI/I) associated with rotation (since 
AI is proportional to the square of the gyro's angular velocity), which leads to a time 
variation in the polhode periods. In addition, we note that the two different quadrupole 
moments also affect the orbit equations resulting in a concomitant contribution to the 
spin precession rate [20| . The inter-connection of these various parameters presents a 
challenge to the interpretation of the data, particularly for the determination of fi^s, 
which is 161 times smaller than ilgs- The data analysis is further aggravated [HI [55] 
by the fact that unexpected classical torques are produced on the gyros, together with 
damping of their polhode motion, which are caused by interactions between the gyro 
rotors and their housings due to electrostatic patches (potential differences). 

Since nearly five decades has passed since the GP-B experiment was conceived and 
designed, there has been much progress in achieving the same results by a variety of other 
methods which we will discuss below. In particular, there are now three independent 
determinations of Qgo consistent with the predicted results arising from Einstein's theory 
and, in fact, one of them pertains to a 2-body system for which additional interesting 
features arise. Also, whereas it is reported that GP-B has confirmed the predicted spin- 
orbit result to better than 1% [H], Murphy et al. [TB] are claiming 10 times greater 
accuracy based on lunar lasing. Thus, more attention will be focussed on the main 
"selling-point" of the experiment, that is the determination of Qgs- But, whether this is 
achievable is very problematic in view of the polhode problems already mentioned and 
other various classical torques on the gyroscopes [22] . Of course, there have been many 
significant "spin-offs" connected with the experiment. But, whether or not this will be 
regarded as justification for the time spent and the cost (over $760 million) and in view 
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of the fact that it required the U.S. Congress to keep it ahve after NASA "-has tried to 
cancel the mission at least three times-" [Ml ES] will surely give rise to an analysis of 
the sociology and politics of science research and funding, along the lines that Collins 
presented for the search for gravitational waves [5S] . 

3'2. Determination of flso from lunar-laser-ranging measurements. - As already al- 
luded to above and discussed in detail in [TO], the de Sitter precession Qgoi^ds) is also a 
frame-dragging effect and "-provides an accurate benchmark measurement of spin-orbit 
effects which GPB needs to emulate-." [T^. In essence, starting with the work of de 
Sitter and investigated in detail by many groups [271 [151 [16] , it is clear that "-because 
of its distance from the sun, the earth-moon system can be regarded as a single body 
which is rotating in the gravitational field of the sun. In other words, the earth-moon 
system is essentially a gyroscope in the field of the sun and its frame-dragging effect due 
to interaction with the sun has been measured, using lunar laser ranging-" [12]. First 
we note that the relevant Ug is 9.8724 x 10~^ and w = 1.991 x 10~'^s~^. The product 
^ttgiD results in a de Sitter [HI [12] spin-orbit contribution which is reflected in a con- 
tribution to the lunar perigee ~ 19.2 m-sec/yr [2T. Second, and more recently. Murphy 
et al. [16| pointed out that the current accuracy of lunar lasing ranging is such that it 
provides a test of spin-orbit coupling to "-approximately 0.1% accuracy - better than 
the anticipated accuracy of the gravity - Probe B result-" and that "-a new effort in 
LLR is poised to deliver order-of-magnitude improvements in range precision" . 

3'3. Ciufolini's Lageos experiment. - The fundamental idea for this work goes back 
to Ciufolini's work in mid 1980, details of which may be found in [28]. In their recent 
paper [TUlllS], Ciufolini and Pavlis used in their analysis two satellites, Lageos (launched 
in 1976) and Lageos2 (launched in 1992) which simply consist of a heavy sphere of retro- 
reflectors which reflect short laser pulses sent from earth, leading to a determination of 
earth-satellite distances with a precision of a few mm. The focus is on the orbital motion 
of the satellites and speciflcally a determination of <Kl/dt, the nodal precession, due to 
the spin of the earth. As already emphasized, this is a spin-orbit effect, involving the 
spin of the earth and the angular momentum of the orbit; the only difference with the 
corresponding spin-orbit measurement in the GP-B experiment is that the spin under 
discussion there is the spin of the gyroscope (the spin of the earth coming into play only in 
the spin-spin interaction involving the earth spin and the gyroscope spin) . Our emphasis 
on precise terminology stems from the fact that there is confusion associated with the use 
of the terminology " Lense-Thirring effect"; some authors, such as [10], use it when the 
earth's spin is playing a role whereas others, such as [3 [22] use it to denote the interaction 
between the earth's spin and the spin of the gyroscope. The emphasis in [TU] on nodal 
motion is based on the fact that, in this case, non-gravitational perturbations are easiest 
to handle. Also, the use of two satellites made it possible to eliminate the effect of the 
earth's quadrupole moment. The use of the GRACE earth gravity model [29], resulting 
in very small uncertainties arising from other harmonics, led to a value of AS>.2masyr~^^ 
for the nodal precession, accurate to a precision of 10%. This is a very impressive result 
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and the error estimate is bound to become even smaller as better results emerge for the 
various multipoles of the earth's field and for the expected improvements in laser-ranging 
precision. For an excellent review of recent measurements of frame dragging using earth- 
orbiting satellites, we refer to Ciufolini's Nature review [33], where he remarks that the 
laser-ranged Italian satellite LARES should " - - in future provide an improved test of 
the Earth's gravitomagnetism with accuracy of the order of 1%". 

4. Two-Body Systems 

4'1. The binary pulsar PSR 1913+16 Q'. - It was the first discovery of a pulsar in a 
binary system and stimulated Barker and me to extend our one-body analysis |3| to the 
two-body arena [8] , the results of which were immediately applied to the Hulse- Taylor 
system [S]. 

In the two-body system, new and unexpected results were obtained. The only previous 
two-body calculation was that of Robertson [3T] who neglected spin and considered only 
two test particles for which he calculated the periastion precession, the result for which 
was the same as the one-body calculation except for the replacement m2 — > ?ni+7Ti2 = M . 
However, when spin is taken into account, the situation is very different. What we found 
is that, for the spin precession of body 1, the replacement in the formula for ^iJ is 
7712 -^ 1TI2 + (/^/3) and, in the case of illo , the replacement is rrii -^ nii + (/i/3), where we 
recall, from (|2.ip . that /i is the reduced mass. Similar replacements occur for the orbital 
precessions H,*^^^ and D,*^'^\ associated with the spin-orbit interaction contributions from 
the spins of bodies 1 and 2, respectively. This is shown explicitly in (|2.17p to (|2.19p . 
Here Ofg = 2 x 10^^, compared to a value of 7 x 10^^*^ for the earth-gyro system and thus 
we expect correspondingly large numbers for the various precession angles. In line with 
these expectations, the observations give il*^^") = 4".2/yr, in agreement with theory. 
The theoretical value for Qlo is 1^.1/yr (where nii = 1. 42771© refers to the spinning 
pulsar) but this quantity is difficult to deduce from the data. Two obvious ways suggest 
themselves: (a) wait for the pulsar beam (which is directed in a narrow cone toward the 
earth) to go outside the line of sight, (which could take a long time since the precessional 
period is about 327 years), or (b) observe the changes over time in the shape of the 
pulses, which have now been seen [321 I33j and which give qualitative agreement with 
theory but due to uncertainties in the beam shape, a quantitative comparison has not 
yet been possible. 

4'2. The binary pulsar PSR B1534+12 flT^ - - Here, based on the two-body spin 
precession formula given in [H] and [5], the theoretical precession rate is O.Sl^/yr. Obser- 
vations amounting to over 400 hours led to a determination of the time evolution of (a) 
the angle between the spin and the magnetic pole and (b) the minimum impact angle of 
the magnetic pole on the line of sight, from which the authors deduced direct evidence 
of geodetic (spin-orbit) precession. This was bolstered by an independent determina- 
tion based on time variation of the shape of the profile of the pulses. The conclusion 
yielded a "-measurement of the precession time consistent with the predictions of general 
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relativity-" fTf and the authors note that "-ahhough the precision is, as yet, Hmited-", 
future plans should allow for improvement. 

4-3. The double binary system PSR J0737-3039 f23 \3^- - This is a highly significant 
discovery, due to the fact that it is the only binary system discovered in which both 
neutron stars are radio pulses; in particular it widens the scope for tests of general 
relativity 35] and for the study of the magnetospheres surrounding pulsars. This binary 
has an orbital period of 2.4 hours and the periastion precession was observed to be 
16°.9/yr. Pulsar A has a period of 22.7ms whereas pulsar B has a period of 2.77s, that 
is A spins 122 times faster. A detailed listing of measured and derived parameters are 
listed in [181 [34l [35] ; in particular, in units of a solar mass, M = 2.58708 with estimated 



values of ruA = 1.3381 and ms = 1.2489. This leads to ag w 4.4 x 10~^. In addition, 
the corresponding reduced mass /i is 0.64596. Thus, using the two-body spin precession 
formula [Hl[l], given in (|2.17p and (|2.18p . the excepted spin-orbit spin precession rates are 
4°.8/yr for pulsar A and 5°.l/yr for pulsar B (which corresponds to geodetic precession 
periods of 75 years for A and 71 years for B). Since the corresponding precession rate 
for the GP-B earth is 1.838 x 10^'^ degrees per year, we see that the number for pulsar 
A is larger by a factor of 2.7 x 10'^. Despite these large values, they have not yet been 
observationally verified. However, it may be possible to deduce the effect of the two-body 
spin-orbit term from the fact that it also results in a value of « 4.06 arc-sec/yr for the 
rate of precession of both the angular momentum and the Runge-Lenz vector of the orbit 
about the pulsar spin direction 36J . This is reflected in a precession of the inclination 
of the orbit [36j and there is also a contribution to the perihelion precession. However, 
in the former case, such an observation is not helped by the fact that the orbital plane 
is close to the line of sight, the angle of inclination being 88°. 69 [35]. The latter case 
is thus probably more promising since the perihelion precession has been measured to 
a relative precision « 10~^. We expect that the second-order spin-independent post- 
Newtonian contribution to the periastron precession to be smaller than the lowest-order 
Einstein- Robertson contribution by a factor « ctg = 4 x 10~^. On the other hand, the 
lowest-order spin-orbit contribution of the fast pulsar A will be smaller by a factor, F 
say, where 



(4.1) ..^(5^^ 

where uj^^' — 44.0545s^^ is the precession frequency of pulsar A [35j and where R w 
15kms is the radius of the pulsar A neutron star. Thus, since {uj^^^ /Cii) = 3.89 x 10'^, 
and (R/a) ~ 1.71 x lO"'^, we see that i^ = 1.1 x 10^*. Thus, the lowest-order spin-orbit 
contribution could be « 10 times the second-order spin-independent contribution and 
also the relative precision « 10^^. This translates to a contribution to the periastron 
precession w (16.9 x 10"'*) jyr k, 5 arc-sec/yr, keeping in mind the uncertainties under- 
lying this calculation. However, if the observations provide even better relative precision 
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in the future, they could "-provide a way to obtain accurate information of the moment 
of inertia of neutron stars-." [36] . 

5. Conclusions 

On the theoretical side, the general relativistic Hamiltonian for a two-body system 
with spin, to order c~^, contains spin-orbit and spin-spin contributions, analogous to 
the corresponding Hamiltonian for QED. These terms separately manifest themselves as 
spin and orbital precessions, which are in fact intimately related, reflecting their common 
origin. Also, it was convenient from an observational point of view, to separate the orbital 
precession into periastron, nodal and inclination angle precessions. Thus, a measurement 
of any of these precessions is potentially a test of gravitomagnetism (spin-orbit and spin- 
spin effects). 

As we discussed above, there is already strong evidence (subsects. 3.2, 3.3 and 4.2) 
for the predicted spin-orbit effects in general relativity. This bears out the conclusion of 
a NASA 2003 scientific review who "-found "some erosion" in the scientific value of the 
frame-dragging experiment-" |37| . However, it does appear, from the initial announce- 
ment of specific results from the GP-B experiment that an accuracy better than 1% may 
have been achieved for the geodetic effect but the lunar lasing experiment is already 
much better, achieving an accuracy of 0.1%. Also, the other experiments, as distinct 
from GP-B, have the advantage of being on-going experiments. For readers interested 
in other theories of gravitation (for which I believe there is neither strong theoretical, 
experimental or observational evidence), a detailed discussion is given in [38j . In par- 
ticular, it is clear that in all these other theories, the intimate relationship between the 
spin-orbit measurements of the GPB group and those of Ciufolini and Pavlis [T^ still 
holds. 

We emphasize that our considerations were confined to order c~^, which is an excellent 
approximation for the systems considered for which the largest value of ag is 4.4 x 10~^. 
However, motivated by the search for gravitational radiation, there is widespread activity 
in investigating the final evolution of a binary spinning black-hole system; the results 
given in 8J provide a touchstone for sufficiently large separations but strong field effects 
soon play a dominant role as the merger takes place so that massive computational effects 
are necessary to extract meaningful results [35] • 

Finally, we note that the theoretical work underpinning the results presented here are 
treated in depth in a companion paper [40] . In particular, we discuss relevant conceptual 
matters dealing with the relation between velocity and momentum (especially for parti- 
cles with spin), the non-uniqueness of spin supplementary conditions and the choice of 
coordinates even at the classical special relativistic level, the fact that a spinning particle 
necessarily has a minimum radius, the corresponding concepts in quantum theory (relat- 
ing to such topics as the Foldy-Wouthuysen transformation and position operators) and 
the fact that the spin effects in quantum electrodynamics (obtained from one-photon ex- 
change) have their analogy in general relativity (obtained from one-graviton exchange or 
purely classical calculations) to the extent that except for (important) numerical factors. 
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the latter results may be obtained from the former by simply letting e^ —^ Gmim2. 

6. Post-Script 

A recent Comment by Kopeiken [41] argues that " - - lunar laser ranging (LLR) is 

not currently capable to detect gravitomagnetic effects " , based on the fact that a 

Newtonian-like translation from the solar system barycentric (SSB) frame to the geocen- 
tric frame (that is the frame co-moving with the earth) gives rise to additional gravit- 
omagetic terms which " - - makes LLR insensitive to the gravitomagnetic interaction" . 
As an example of such extra contributions, we note that coordinate tranformations give 
rise to Thomas precession contributions of the same order of magnitude as those arising 
from post-Newtonian effects [131 [55]. However, for the 3-body earth-moon-sun system, 
there are several velocity-dependent contributions and a Reply by Murphy by et. [33] to 

the Kopeiken Comment disputes his conclusion on the basis that " transformations 

of the velocity-dependent terms from one frame to another [are] strongly constrained 

by experiment". 

More recently, Ciufolini in a lucid analysis [33] , reiterates a point which he made earlier 

in his book with Wheeler [2H] that " the de Sitter effect and the Lense-Thirring drag 

are fundamentally different phenomena " |28^ since the curvature invariant *R.R is 

zero in the former case (which is generated by the translational motion of a mass) but 
non-zero in the latter case (which is generated by the rotation of a mass) . Ciufolini then 

goes on to stress that " one cannot derive rotational gravitomagnetic effects from 

translational ones, unless making additional theoretical hypotheses such as the linear 
superposition of the gravitomagnetic effects." While we agree with Ciufolini, from our 
point of view, both of the effects under discussion arise from a common spin-orbit term 
in the basic Hamiltonian H, which can manifest itself in various ways. In particular, as 
detailed above, the same spin-orbit term gives rise to both a precession of the spin (or 
rotation axis) per se and also a precession of the perihelion of the orbit. In particular, 
the rate of precession of the spin arising from the spin-orbit term in H depends on the 
angular momentum and not on the spin whereas the rate of orbital precession depends 
on the spin but not on the angular momentum. The rate of precession of the GP-B 
gyroscope orbiting the earth is representative of the former case and so is the earth- 
moon gyroscope orbiting the sun. Concomitantly, in both cases, there is an associated 
orbital precession; this is negligibly small in the case of the GP-B gyro but not so for the 
case of the earth-moon gyroscope because the latter is itself a two-body system forming 
part of a three-body system. As a consequence, the motions of the earth and the moon 
can be separately tracked leading, in particular, to spin-orbit amplitude contributions to 
the lunar orbit at the 6 m level [IB] . 

We should also emphasize that, as already discussed, there are two separate spin-orbit 
contributions. However, in all one-body problems, the dominant spin-orbit contribution 
to the spin precession is associated with the spin of the smaller mass [see (2.17) and 
(2.18)] i.e. the gyro mass in the case of GP-B and the earth-moon mass in the case 
of LLR. Of course, there is also a spin-spin contribution in both cases; an attempt 
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has been made to measure this in the case of GP-B but it is generaUy negUgible for 
contributions to one-body orbital precessions (compared to spin-orbit contributions and 
the dominant Einstein- Robertson periastron precession term). On the other hand, for 
two-body systems, such as the double binary system, both spin-orbit contributions to 
spin precession are important and, as already noted, in Section 4.3, their contribution to 
the periastron precession of the orbit of the double binary system may also be measurable 
in further observations. 

More generally, we emphasize that our analyses and calculations demonstrate the 
universality of gravitational spin effects regardless of whether the source is the spin of an 
elementary particle, the rotation of a macroscopic body or a translational effect. I thank 
Professor Herbert Pfister for comments which initiated this post-script. As a further 
update, we note that our result for spin-orbit precession in a 2-body system [S] has been 
verified, to an accuracy of 13%, by Breton et al., who observed a precession of the spin 
axis of pulsar B in the double binary system by an amount 4.77°/yr. [45j . 
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